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Abstract. We perform a Bayesian model comparison for scenarios within the qnadratic 
cnrvaton model, determining the degree to which both are disfavonred with respect to the 
ACDM concordance model and single-field qnadratic inflation, nsing the recent Planck data 
release. Despite having three additional model parameters, the simplest cnrvaton scenario is 
not disfavonred relative to single-field qnadratic inflation, and it becomes favonred against 
this single-field model when we inclnde the joint BICEP/Keck /Planck analysis. In all cases we 
assnme an instantaneons inflaton decay and no snrviving isocnrvatnre pertnrbations. Despite 
the snccess of Planck reaching its forecast measnrement accnracy, we show that the cnrrent 
constraints on local non-Ganssianity are insnfhciently precise to have any significant impact 
on the evidence ratios so far. We also determine the precision ct(/nl) reqnired by fntnre 
measnrements assnming a fidncial valne of /nl = —5/4 or 10.8 to no longer disfavonr the 
cnrvaton against the ACDM parametrisation, and we discnss the effect that the predicted 
increase in precision from fntnre measnrements on /nl may have. We show that onr resnlts 
are not very sensitive to onr choice of priors. 
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1 Introduction 

Following the recent Planck data release, which has not detected any clear deviation from the 
predictions of the simplest single-field inflationary models, it is interesting to ask whether this 
should be taken as evidence that inflation was driven by a single field, or whether the data 
is still not good enough to discriminate between single and multiple-field models of inflation. 
In general, more complex models may have a better fit because they have more parameters 
which can be tuned, but they should be penalised if large amounts of parameter space are a 
bad fit to the data. 

We study the curvaton scenario, which is one of the most popular and widely studied 
multifield models of inflation. In the original scenario [1-3] (see also [4, 5]), the curvaton field 
generates the primordial curvature perturbation while the inflaton field drives inflation. By 
definition, the curvaton field is light and subdominant during inflation, but its energy density 
grows relative to the background energy density after inflation ends, due to the curvaton 
decaying much later than the inflaton. The curvaton scenario is perhaps the best known sce¬ 
nario for generating large non-Gaussianity (of the local shape), a potential signature which 
Planck has constrained much more tightly than the WMAP satellite but not detected. Build¬ 
ing on the original scenario, one can also consider mixed scenarios in which both the inflaton 
and curvaton significantly contribute to the power spectrum [6, 7], or the inflating curvaton 
scenario [8-11], which changes the observational predictions of the model. 

We perform the first Bayesian analysis of the “simplest” curvaton scenario, consisting of 
an inflaton and a curvaton field both with quadratic potentials. This model was originally 
studied in 2002 [12], an analysis of the observables produced was performed in 2008 [10] and 
the predictions of the model were revisited following the BICEP2 data release in 2014 [13]. 
While those papers studied the predictions and parameter space of the model, there has been 
no work asking whether the data prefers the addition of the curvaton field. Interestingly, 
although this model is famous as a means of generating large non-Gaussianity, we find it is 
not the non-detection of non-Gaussianity which puts this model under observational pressure. 

Recent papers have conducted comprehensive model selection analyses over parameter 
spaces from single-field inflationary potentials [14-18] yet when reviewing the literature, there 
appears to be very little using the Bayesian paradigm over multi-field models such as the 
curvaton. We note that a Bayesian analysis with WMAP data of a Peccei-Quinn field has 
previously been performed in 2004 [19] which is effectively a physically motivated mixed 
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inflaton-curvaton generation of the curvature perturbation. Best-fit parameter values over 
the {ns,r} space have also been previously obtained in [20] and sampling over the same space 
in [9], both for different choices of prior. 

Conducting a model comparison on the curvaton presents a interesting test of how 
discerning the current data is, given that the choice of prior on parameters within the model 
can be difficult and misleading in certain cases. A particular choice of prior to carefully 
consider will be that of the vacuum expectation value of the curvaton cr*, which we choose to 
have either a Gaussian or flat prior — a physically motivated or statistically relevant choice 
respectively — and it is important to discuss how independent of these prior choices our 
conclusions are. 

In this work; we make the initial step in discerning how a model selection on the curvaton 
is to be accomplished; we perform a Bayesian analysis to ascertain how much the data may 
prefer, if at all, the single-field quadratic inflaton over adding a quadratic curvaton; and we 
quantify how useful future measurements on local-type primordial non-Gaussianity [21-24] 
will be to favour or disfavour the curvaton. 


2 The curvaton scenario 


We shall firstly review the general predictions of the curvaton model with potential [12] 
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that can be used to make predictions of 4 specific observational parameters: the power spec¬ 
trum amplitude "P^, the scalar spectral index n^, the tensor-to-scalar ratio r and the non¬ 
linearity parameter /nl which characterises the degree of local type non-Gaussianity. 

To find the total power spectrum, one simply sums the inflaton and curvaton 
contributions 


= Vf + V? = 
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where hereafter denotes the value as the current scale crossed the horizon during inflation, 
Mp is the reduced Planck mass and r^ec is defined as [1] 
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where ‘dec’ denotes evaluation at curvaton decay time. In the original curvaton model one 
assumes that is negligible, we here are more general and allow an arbitrary contribution 
to the power spectrum from both fields, significantly increasing the available parameter space 
for curvaton-like scenarios [6, 7]. 

The slow-roll parameters are 
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where and is the number of efolds during inflation from when the current scale 

crossed the horizon until the end of inflation (denoted by a subscript ‘end’ below), which is 
given by 
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Note that the curvaton provides a negligible contribution to e because the inflaton dominates 
the total energy density as it drives inflation. 

Using the slow-roll parameters one can hnd the spectral index [25] 


n, - 1 = ^(-2e + 2r?,) + ^(-6e + 2r ]^), (2.6) 

derived by parameterising the degree to which the inflaton fluctuations contribute to the total 
(i.e. observed) power spectrum amplitude V(^. For a given value of e, the tensor-to-scalar ratio 
can only be suppressed compared to its single-field value, and is given by [25] 


r = 



(2.7) 


Local non-Gaussianity is a characteristic observable of the multiheld inflation models 
which if detected would rule out all single-held inhationary models [26] (however see [27-29]). 
The curvaton model can produce 10^ > /nl > —5/4 partially depending on how late the 
curvaton decays, parametrised by rjeo where [13] 
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where rjec = 1 acts as an attractor for the majority of points within its parameter space, 
because rjec initially grows quickly and then asymptotes to unity once it dominates the 
background energy density, unless the curvaton decays early enough. We therefore hnd /nl = 
—5/4 in the limit of <C V(^ and a sufficiently late decaying curvaton. For the choice of 
priors we make (which are broad), we hnd the majority of non-Gaussianity values are in the 
region 10 > /nl > —5/4 and are hence not in tension with observations. 

Implicit in all of the expressions above, W has a general expression in all quadratic 
curvaton models. By assuming instantaneous inhaton decay [13] 


1 1 

W = 58 + -ln^--iV^at, (2-9) 

4 4 

where W is taken to correspond to the horizon crossing time of the Planck pivot scale, 
/c* = 0.05Mpc“^. The expression includes A^mat) the number of efolds with a matter-like 
equation of state caused by the curvaton held, while it oscillates after becoming dominant. If 
the curvaton never dominates the background energy density then iVmat = 0. 

The key value of r^ec = 3/7 corresponds to inhaton-curvaton equality, p^f, = pa at the 
curvaton decay time, at which time the Hubble parameter equals the curvaton decay rate. 
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To- = H\dec- We find two regimes are possible. Firstly, for Fq- > F^eq. the curvaton decays 
before it can dominate the total energy density and = 0. Secondly, using the relation [12] 


^dec 





( 2 . 10 ) 


we can obtain an expression for At^at in the opposite regime Fg- < FFeq.) which means that 
the curvaton decays after it dominates the energy density with r^ec >3/7 and 
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We note that by setting » in all of the expressions above one retrieves the 

predictions of the single-field quadratic inflaton which becomes a limiting case of the curvaton 
model in which no significant perturbations are generated by the curvaton but may still 
vary through equation (2.9) and A^matj which is zero without the curvaton. The curvaton 
models are therefore found by setting arbitrary power spectra of the inflaton and curvaton in 
the equations above, yielding 3 interesting scenarios in particular: 

1. The mixed inflaton-curvaton scenario allows free choice of any sensible M so that 
both inflaton and curvaton fluctuations may significantly contribute to the power spec¬ 
trum. This introduces a free effective parameter 

2. The pure curvaton scenario is a subset of the mixed inflaton-curvaton scenario, match¬ 
ing the original model proposed by [1]. The curvaton is the sole generator of the fluc¬ 
tuations, i.e. Vq <C Vq. 

3. The dominant curvaton decay scenario is a subset of the pure curvaton scenario 
which fixes rjec = 1 (/nl = ~5/4) where the curvaton decays only after it has become 
dominant in the energy density, rjec = 1 is the most likely value for both of the curvaton 
scenarios considered above, so this is an interesting case to consider. 

It is important to note that in all of the scenarios above we assume that there are no 
surviving isocurvature modes, which in practise means that either the curvaton dominates 
before decay and/or that the Universe thermalises after the curvaton decay. Whether the 
isocurvature perturbations survive until the CMB forms depends on the microphysics of re¬ 
heating and this is beyond the scope of our investigation. At any rate, the assumption that 
the Universe fully thermalises necessarily ensures full isocurvature decay. 

We note that throughout this paper we make the common simplifying assumption that 
the inflaton does not decay into the curvaton after inflation. However several papers have 
shown that even if only a small fraction of the inflaton field decays into the curvaton [5, 30] 
then the amplitude of the curvaton perturbations is suppressed and that the value 

of /nl may be enhanced. Including this effect goes beyond the scope of this paper, but the 
increased freedom in /nl may make the Planck non-Gaussianity constraints more effective in 
this case. 

Fig. 1 illustrates the extrema boundaries available to the curvaton scenario over {ns,r} 
space in this paper. The single field inflaton is represented by the larger black point cor¬ 
responding to Nj, = 58 from equation (2.9) which corresponds to a single point because we 
assume instantaneous decay of the inflaton. As shown in Fig. 1, the tensor-to-scalar ratio 
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can be larger than the single-field valne, in the case where the cnrvaton pertnrbations are 
negligible bnt its effect on the backgronnd evolntion generates a large A^mat, which can rednce 
to a minimnm of A^* =48, making e larger. The mixed inffaton-cnrvaton scenario is fonnd 
from all points within the region defined between the three black lines. The predictions of 
the pnre and dominant cnrvaton scenarios lie along the bottom of the plot (r ~ 0) between 
the two black lines. 

We note that Fig. 1 is derived nsing the single field consistency relation for r ~ this is nsed 
to evalnate r with a pivot scale = 0.002Mpc“^ - bnt by contrast the scalar spectral index 
Us has no associated scale becanse the rnnning is assnmed to be 0. This leaves ns with a small 
fnndamental error in the exact placement of the point over {ns,r} for the single field model 
becanse it predicts that both Ug and r change eqnally with the pivot scale (n^ oc r oc 1/A^*). 
The cnrrent constraint on the scalar rnnning dng/dlnk = —0.003 ± 0.007 [14] is close to 0, 
therefore we do not expect this effect to be significant to onr overall analysis. In Fig. 1 we 
calcnlate Ug and r at the nsnal Planck pivot scale of O.OSMpc”^ nsing (2.9), which is the best 
choice for the spectral index bnt does introdnce a small error into r when it is large. 

3 Model Selection 

3.1 Selection methods and prior ranges 

In the attempt to qnantify the degree to which the cnrrent data is able to discern between 
models, one might expect the cnrvaton to be largely disfavonred against a single-field potential 
nsing the Bayesian approach. The cnrvaton boasts two additional parameters to its potential, 
pins the late time parameter rjec with no fnrther constraints, thereby losing ‘efficiency’ in 
statistical fitting. 

We shall show that the cnrvaton performs mnch better than one might expect, highlight¬ 
ing some of the weakness in cnrrent data so far in discerning between models. In particnlar, 
the cnrvaton is not disfavonred against single-field quadratic inflation. To this end, we now 
briefly review the model selection methods nsed here. 

We maximise the model likelihood £ taking into acconnt the data provided by Planck [31] 
to find ^xlff 

Axlff = 2[ln£max(Mi) - ln£max(A42)] . (3.1) 

Models which fare well in a comparison nsing ^Xeff be as statistically favonrable 

in the Bayesian paradigm dne to a wastage of parameter space penalised by the evidence £ 
with model parameters k 

T = y d’^xPT{x\M.)C{data\x), (3.2) 

in ratio with another reference model and is typically qnoted logarithmically. Once calcnlated, 
the ln{£/£ref) valnes mnst have some measnre to correctly interpret their relevance in a model 
comparison, one snch scale nsed in [32] is that of Jeffrey’s which is shown in Table 1. We are 
carefnl here to note that the Jeffrey’s scale cannot be nsed as mnch more than a gnideline, 
especially in instances where models appear to have similar evidences as it can be shown to 
be an incomplete comparison tool [33]. 

An important concept with the introdnction of Bayesian statistics is the choice of pri¬ 
ors, which may be categorised as either physically motivated or non-informative (Jeffrey’s 
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Figure 1. The {ns,r} plot taken from [14] outlining the regions where the larger dot is the inflaton 
single field inflaton limit at fV* = 58, the black lines contain the mixed inflaton-curvaton and the r = 0 
values within this region belong to the pure curvaton scenarios. The left black line that corresponds 
to iV* = 48 at the top is essentially straight due to an approximate cancellation between the second 
and third terms of equation (2.9). The right hand black line follows by maximising iV» and m, whose 
upper limit we set as m = M/2 so that the curvaton does not roll significantly during inflation and in 
agreement with [13]. However the plot differs from the equivalent one in [13] because we here calculate 
the allowed range of iV* as a function of all the model parameters, rather than assuming a fixed range 
of Njf. The Planck data used in this paper corresponds to the blue contours and region. 


\H£/£ref)\ 

Interpretation 

< 1 

Inconclusive 

1 - 2.5 

Weak evidence 

2.5 - 5 

Moderate evidence 

> 5 

Strong evidence 


Table 1. The adapted Jeffrey’s scale as used in [32] that assigns meaning to \n{£/£ref) values. 

prior) [34], the latter of which is proportional to the square root of the Fisher information 
determinant. 
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Model Scenarios 

Prior ranges and constraints {Mp = 1) 

Single field inflaton 

—4 > log;^o("^) > —6.75 [15] 

Mixed inflaton-curvaton 

-4 > logio(M) > -15.3 [15] 
logio(M/2) > logio(m) > -39 [12] 
(flat or Gaussian) a* < 0.01 [7] 
logio("i) > logio(r^) > -39 [12] 

Pure cnrvaton 

Mixed inflaton-curvaton subset with extra constraint: 

Vf < 0.01 (P^ + V^) = O.OlPf 

Dominant cnrvaton decay 

Pure curvaton subset with extra constraint: 

^dec 1 

^ logio(m) -b 41ogio(o-*) - 4 > logio(r^) (2.10) 

ACDM concordance 

1.02 >ns> 0.9 [15] 
3.2 > In(lOiop^) > 3.0 [15] 
r = 0 [15] 


Table 2. The list of prior ranges for all of the inflationary scenarios discnssed in the text. In 
all models we set an overall prior of 10“^ > Vt^ > 10“^^, emnlating [15]. Unlike [15] we do not 
set a prior on Afj, becanse it can be precisely calcnlated after all of the above model parameters 
have been specified. In the case of the cnrvaton scenarios we have placed an additional constraint of 

10 -'^ >vi> 10 -“. 


We present our prior ranges for the scenarios discussed previously with our rationale in 
Table 2, where we shall now discuss our choices. The single field quadratic model in Table 2 
sets a non-informative logarithmic prior only on M, matching that of [15] as we find that 
W = 58 from equation (2.9) for an instantaneous decay scenario (which means that the 
Hubble crossing value of 4> is not a free parameter). 

In all cnrvaton models we set a logarithmic prior over M to match the single-field 
inflation model. In the mixed inflaton-curvaton scenario, M is given a full range between 
its upper limit corresponding to the single field inflaton and its lower limit to the minimum 
M allowable by the prior 10“^ > V(^ > 10“^^ which we attribute globally to all models 
keeping only the region of parameter space that is compatible with basic structure formation 
requirements [15]. In contrast, the pure cnrvaton must have a tighter upper limit on M (by 
a factor 10 compared to the mixed and single-field upper limits) resulting in a suppression 
of r. Because the choice of prior is logarithmic over the inflaton mass we find the effective 
parameter volume available to the pure cnrvaton scenario to be not much smaller than the 
mixed case. 

In all relevant cases we also use a non-informative logarithmic prior over m and obtain the 
upper limit from the requirement that the cnrvaton does not roll significantly until inflation 
ends, requiring m < M/2 which implies r]„ < rj^/A. In choosing the lower limit on m and Tq- 
one can allow for variations of baryogenesis to offer constraints however we choose the most 
conservative bound, which is big bang nucleosynthesis where Hnud. — 

We set a firm upper limit of |(t*| < IQ-^Mp to stay within the regime where our equations 
for the cnrvaton are valid [7] and to avoid the inflating cnrvaton scenario. We have two choices 
of prior on cr*. As the order of magnitude of u* is much better known, the non-informative 
prior one can place is a flat prior. We also have a physically motivated alternative following [35] 
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to the obtain the equilibrium distribution of cr* in a de Sitter background 

Pr((T*) oc exp 

This distribution is obtained due to the relaxation of the vacuum expectation value at the equi¬ 
librium limit being normally distributed over fi* = 0 with a standard deviation oc M'^N^/m. 
We therefore carry out our analysis separately with both flat and Gaussian priors over cr* 
to assess the effect, if any, this may have on the evidences. Note that even when using a 
Gaussian prior on cr*, we truncate the range of the Gaussian to enforce |(T*| < 

Finally, the prior we place over Fq- in all scenarios is logarithmic where we hnd the ranges 
differ, setting an upper limit for the mixed and pure cases as m so that the curvaton oscillates 
before it decays and an upper limit for the dominant curvaton scenario from hxing rjec = 1 in 
equation (2.10). The lower limit of in all 3 curvaton scenarios follows from the constraint 
that the curvaton must decay before big bang nucleosynthesis. 

3.2 Calculating the evidences 

The evidence is far more computationally expensive to calculate than the maximum like¬ 
lihood, particularly with large numbers of parameters, so other papers [14-17] have imple¬ 
mented nested sampling techniques to compensate for the increase in parameter volume. By 
comparison direct Monte Garlo methods - sampling across a grid of sufficiently narrow step 
size - are still feasible as a way of estimating the evidence for relatively small parameter 
spaces such as {Vc^, Ug, r, /nl} and we shall implement these methods in this paper. 

Initially, we choose the following computation ranges and bin widths: 

3.4 > ln(10i°P^) > 2.8, 1.02 >ns> 0.9, 0.5>r>0, 30 > /nl >-20, 

Aln(10^°P^) ~ 0.006, 0.001, Ar~ 0.005, 

where by marginalising only over the /nl measurement we reduce the effective number of 
parameters to be stored so that the data distribution can be represented in a multidimensional 
array and a Gaussian function is placed over /nl = 0.8 with cr(/NL) = 5.0 [36], therefore it 
does not have an associated bin width. The widths are roughly consistent with selecting 
histogram bin sizes following the Freedman-Diaconis rule [37], which takes into account the 
sample size of the Planck Markov chains ~ 2 x 10^ obtained from [31]. 

Implemented for all models in this paper, our numerical method was as follows: 

1. Smooth the binned MGMG chains over {ln(10^*^P(^),n^,r} by kernel. 

2. Multiply the smoothed likelihood by a Gaussian /nl function to obtain the likelihood 
over {ln(10^°P^),ns,r,/NL}- 

3. Obtain 10® sampled points from om code to emulate a prior, where the single held was 
an effective delta function. 

4. Sum the resultant likelihood from 2. over the prior mass to numerically approximate 
the evidence integral. 
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Model Scenarios 


\ci{E/Eref) 

Single held inhaton 

6.6 

-5.2 

Mixed inhaton-curvaton (cr* hat prior) 

3.7 

-4.9 

Mixed inhaton-curvaton (cr* Gaussian prior (3.3)) 

3.7 

-4.6 

Pure curvaton (cr* hat prior) 

4.7 

-4.7 

Pure curvaton (cr* Gaussian prior) 

4.7 

-4.6 

Dominant curvaton decay (cr* hat prior) 

4.7 

-4.6 

Dominant curvaton decay (cr* Gaussian prior) 

4.7 

-4.4 


Table 3. The \n{£/Eref) and ^Xeff comparison statistics in the form of ratios against the ACDM 
refererence model over {Vi;,ns,r, f^^,} space with the Planck TT+lowP+BAO dataset as in [14]. 


For a dramatic difference in computational efficiency, when obtaining the smoothed 
likelihood from the MCMC chains, we have chosen to ignore the covariance between the 
background cosmological parameters and our chosen parameters {ln(10^‘^'P^),n^,r} by fore¬ 
going marginalisation. This ignorance necessarily introduces an error into our calculation of 
the likelihood function over {ln(10^*^'P^), r, /nl} and therefore an error in the ln£i ratios. 
To test that this is small, we have compared our evidence ratios of monomial potentials, the 
simple Higgs inflationary potential and radiatively corrected Higgs potential to those obtained 
by Planck 2013, Planck 2015 and [16]. 

The In I? ratios for monomial potentials can be matched to Planck 2013 in [15] up to 
an error of 0(0.1). Planck 2015 use an alternative baseline model {PI? inflation) to calculate 
their evidence ratios, however we are able to verify an error of 0(0.1) when directly taking the 
ratios between our test monomial potentials. Finally, we find our error is also of 0(0.1) when 
comparing our calculated evidence ratios to those obtained in [16] using monomial potentials, 
the simple Higgs inflationary potential and radiatively corrected Higgs potential. 

We present in Table 3 values and evidence ratios computed for each of the scenar¬ 

ios. The results are discussed in the next section. In table 4 we calculate the required 1 — a 
error bar, cr (/nl)) which observations are required to reach for each of the curvaton scenarios 
to have ln{£/Er^f) ~ 0 against ACDM. We do this for two fiducial values of /nl ? 0, firstly 
for the attractor value of (rdec = 1) /nl = —5/4 and secondly for the current 2 — cr upper 
bound value of /nl = 10.8 [36]. These results assume that apart from a detection of /nlj all 
other cosmological data including the observational constraints on the spectral index and r 
would remain the same. 

One can find by (arbitrarily) setting /nl = 0 for all models and re-sampling that the 
ln{E/Ej-ef) values are unchanged up to 0(0.1), which means that despite the curvaton scenario 
being a famous method of generating a potentially large non-Gaussianity, the observational 
constraint on /nl does not significantly disfavour the model or change the Bayesian evidence. 
However, in sharp contrast, by changing instead the power spectrum for all models, fixing 
them to be the same as ACDM concordance with 3.2 > ^(lO^'^'P^) > 3.0, we find an increase 
in the ln{E/Eref) values by approximately 4. 

With the especially tight error bar on V(^ we expect to hnd at least one finely tuned 
parameter in all of the models we study. Such parameters may have wide prior ranges but 
only specific values within those ranges will correspond to a favoured region by the data, 
e.g. logio(Af/Mp) ~ —5.2 for the single held hts the power spectrum well, but the observa¬ 
tional error bar is less than 2% of the full range that is allowable by our globally imposed 
10“^ > V(^> 10“^^ prior. 


9 





Model Scenarios 

hducial /nl 

^^(/nl) 

Mixed inhat on-curvaton (cj* hat prior) 

-5/4 

0.4 

Mixed inhat on-curvaton (cj* Gaussian prior) 

-5/4 

0.4 

Pure curvaton (fi* hat prior) 

-5/4 

0.4 

Pure curvaton (fi* Gaussian prior) 

-5/4 

0.4 

Dominant curvaton decay (a* hat prior) 

-5/4 

0.4 

Dominant curvaton decay (a* Gaussian prior) 

-5/4 

0.4 

Mixed inhat on-curvaton (cj* hat prior) 

10.8 

2.8 

Mixed inhat on-curvaton (cj* Gaussian prior) 

10.8 

2.6 

Pure curvaton (fi* hat prior) 

10.8 

2.4 

Pure curvaton (a* Gaussian prior) 

10.8 

2.6 


Table 4. A list of required /nl measurements with cr(/NL) precisions which observations are required 
to reach for each of the curvaton scenarios to have ln(£/£re/) — 0 against ACDM. The measurement 
at /nl = 10.8 reflects the current 2 — a upper limit from observation [36] and /nl = —5/4 would 
be a strong indication of the curvaton as it is derived from the case where r^ec = 1 and the inflation 
perturbations are subdominant compared to the curvaton perturbations. 


Altering the large ranges of both m and To- by shortening the lower bound from nucle¬ 
osynthesis Mp) to 10“^®Mp had minimal effect on the mixed curvaton log evidence - 

shifting its magnitude by just 0.1 - therefore we hnd with our results against this choice of 
prior is quite robust. 

We set M to its favoured value for the curvaton model in section 3.2 and found essentially 
no change in the evidence, whereas setting the same value for the single held suppressed most 
of the disfavoured region of its parameter volume and therefore the log evidence shifted 
from —5.2 to —1.3. One can explain this difference between models because by imposing a 
logarithmic prior on M, we hnd that most of the parameter volume of the mixed curvaton 
has > 0.99V(^ i.e. the curvaton contributes most to the admixture of the perturbations, 
therefore the parameters that determine its magnitude, rdec and cr* are those most affected 
by hue tuning. 

We further note that because the majority of points in all of the curvaton scenarios 
have rdec = 1) the hnely tuned parameter of the curvaton model is essentially its vacuum 
expectation value fi*. To prove that cr* was hne tuned, we hxed both logiQ{M/Mp) ~ —6.1 
and (T* = 0.01, corresponding to the region in the parameter volume of the mixed curvaton 
that is favoured in the data by and we found a shift in the log evidence from —4.9 to 
— 1.9 as expected. The compensation for the extra parameters in the evidence can therefore 
be traced back to the choice of range in cj*, but is essentially independent of its prior over 
the given range - Gaussian or flat - because the log evidence ratios do not vary signihcantly 
with this choice. 

Performing the same analysis with the chains from the Planck 2013 release [15] we 
found the evidence ratio of the mixed inflaton-curvaton scenario to be more disfavoured at 
—5.2 whereas the single held model was marginally preferred at —5.0. Following on from this 
we also obtained and In^S/S^ef) values of the single held and curvaton scenarios using 

the BICEP/Keck/P/anc^ data from [14] and one can see that the curvaton improves against 
the single-held inhaton, where differences in the log evidence between the two models are 
generally between 1.6 — 2.3. 
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Model Scenarios 

^Xe// 

\n{£/£ref) 

Single held inhaton 

10.4 

-7.0 

Mixed inhaton-curvaton (u* hat prior) 

2.7 

-5.4 

Pure curvaton (cj* hat prior) 

4.8 

-4.7 

Dominant curvaton decay (cr* hat prior) 

4.8 

-4.7 


Table 5. The \\i{£/£ref) and £^x1ff comparison statistics in the form of ratios against the ACDM 
refererence model over r,/ nl} space with the Planck TT+lowP+BKP dataset as in [14]. 


3.3 Discussion 

It is clear from looking at the evidence ratios in Table 3 that the data disfavours both the 
curvaton and quadratic single-held inhation model signihcantly with respect to the reference 
ACDM concordance model where all models listed have close to ‘^strong evidence'” against in 
accordance with Table 1, where the best htting model according to its evidence appears to be 
the dominant curvaton by a fraction that is not signihcant according to the Jeffrey’s scale. 

The data appears to yield no discernible difference in the evidence between the single 
held inhaton and any of the curvaton models, the ^Xeff values indicate that the curvaton 
scenarios ht the data better but are penalised due to their parameter space wastage. The 
mixed inhaton-curvaton scenario has the better ht compared to the other two scenarios from 
its ‘^Xeff value due to the fact that it may cross deeper into the favoured region over {n^, r} 
space with a larger tensor-to-scalar ratio, as can be seen from Fig. 1. 

Performing our analysis with the BICEP/Keck data in section 3.2 hnds a relative im¬ 
provement in the evidence of the curvaton, where on the Jeffrey’s scale, all scenarios are 
categorised to prefer single-held inhation with “ ITea^ evidence”. However, evidence for all 
curvaton scenarios relative to ACDM is almost the same as using Planck data alone, the 
only signihcant change is the tighter constraint on r making the single-held model more 
disfavoured. 

At the end of section 3.2 we calculated the changes that hxing the power spectrum 
predictions to those of ACDM have on the log evidence ratios and it is found in Table 4 
that in all cases the values shift by approximately 4. This is an interesting result as it not 
only demonstrates that the parameter with the largest inhuence on model selection for both 
the curvaton and single held inhaton is but also because the shift is roughly uniform for 
all models, it implies that the difference in the evidences between curvaton and single-held 
inhaton predominately derives from their differing parameter spaces over Ug and r space (but 
not /nl which is too weakly constrained to have much effect). 

The choice of cr* prior appears to have little effect on the evidences of Table 3, the pure 
and mixed curvaton scenarios still are disfavoured versus ACDM concordance with the same 
value and it is only the dominant curvaton scenario that appears to become slightly more 
favoured with the Gaussian prior. The logarithmic priors over a wide range of M and m 
largely reduce the effect that the constraint < O.Ol'Pi;^ has in the pure curvaton scenario, 
making it unsurprising that there is little difference between the evidence ratios of the pure 
and mixed cases. We also found almost no sensitivity to changes in the minimum curvaton 
decay rate, increasing the lower bound on Pq- by ten orders of magnitude led to almost no 
change in the Bayesian evidence ratios. 

In Table 4 we hnd the required measurements and their associated precisions (t(/nl) 
to make the curvaton scenario have the same evidence as the ACDM concordance model. 
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assuming that all other observational data remains the same. We additionally assume that 
the measured value of /nl is centred on the fiducial value, with an error bar given in the right 
column of Table 4. If /nl is as large as currently allowed at 2 — o", then Euclid should reach 
sufficient accuracy to favour the curvaton scenario, but in the limit of the dominant curvaton 
scenario which predicts /nl = —5/4 an error bar significantly below 1 will be required, which 
is forecast to be achievable in the longer term [21-24], Although a much higher precision on 
the measurement of /nl is required if /nl = —5/4, the required significance of the detection 
is about 3-fT, which is lower than the detection significance of 3.5 to 4.5 — u if /nl = 10.8. 
The reason is because a large part of the curvaton prior space corresponds to the dominant 
curvaton scenario, with rjec = 1 and /nl = —5/4. In contrast, only a small region of 
parameter space corresponds to /nl — 10. This also means that if future data does detect 
/nl = —5/4 we would have a strong motivation to believe this comes from the dominant 
curvaton scenario, but if /nl — 10 was detected one would want to also explore other scenarios 
capable of generating non-Gaussianity, to see if any model predicted this value in a large range 
of its parameter space. 

Near the end of section 3.2 we also made a comparison between the impact of the recent 
Planck 2015 release to the previous data from 2013 on the evidence. We find that the mixed 
curvaton becomes marginally more favoured against the single field inflaton by using the 2015, 
shifting from la{£/S^ef) = —5.2 to —4.8. This slight change we deduce is from the tighter 
constraints placed over the tensor-to-scalar ratio in 2015, mildly favouring lower r regions 
than in the 2013 release. 

Notice that by applying the prior |(T*| < 10“^ (a small range compared to the inflaton’s 
initial field value) we do not attempt to quantify whether the curvaton scenario is likely 
within a multi-field context. We are instead asking whether the curvaton scenario is preferred 
within the context of a Lagrangian with two quadratic potentials where one has been given 
a large initial field value to drive inflation, while the second has a small field value and small 
decay rate to potentially act as a curvaton. Within this context, we find that Planck data is 
inconclusive but the addition of BICEP/Keck data prefers the curvaton scenario. 


4 Conclusions 

We find that both single-field quadratic inffation and the curvaton scenario (added to an 
inflaton field with quadratic potential) are strongly disfavoured relative to a ACDM reference 
model, which has Gaussian perturbations and no tensor perturbations. However the data 
is not good enough to discriminate between single-field quadratic inflation and the curvaton 
scenario, despite the curvaton scenario having three additional parameters. The priors are 
listed in Table 2. 

We studied three variants of the curvaton scenario, 1) the original “pure” scenario in 
which the inflation perturbations are neglected, 2) a subset of this case where the curvaton 
dominates the background energy density before it decays and 3) the mixed inflaton-curvaton 
scenario in which the perturbations of the inflaton are allowed to be large. The Bayesian 
evidence was similar for all three cases, probably because the majority of the prior space 
matches the dominant curvaton scenario. The mixed scenario has the best fit, as can be seen 
by eye from Fig. 1, corresponding to a value of r ~ 0.06. However only a small part of the 
available parameter space (were the inflation mass is close to its required value in the single 
field case) generates a non-negligible value of r. In most of the allowed parameter space, r 
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is negligibly small for the mixed scenario, as it is by definition in the pure and dominant 
curvaton scenarios. 

The Bayesian evidence ratios are not very sensitive to our choice of priors. We consider 
both a uniform and a Gaussian prior over the initial curvaton value, fi*, and hnd very little 
change in the results. In both cases we require |cr*| < O.OlMp. Interestingly, we have 
found that the possibility of large non-Gaussianity in the curvaton scenario does not make 
it disfavoured, since if one drops the observational constraint on /nl the evidence ratios 
hardly change. We also apply the priors 10“^ > Vc_ > 10“^^ and 10“^ > > 10“^^. For 

simplicity we assume that the inflation decays instantaneously after inflation ends, dropping 
this requirement would leave the inflation scenario with two free parameters and the curvaton 
with hve parameters. 

The degeneracy in the evidences of single held and the curvaton in this paper demon¬ 
strate a situation where the Jeffrey’s scale has essentially failed to hnd a preferred model 
given the data, despite some penalties included for wastage of parameter space. An impor¬ 
tant next step in order to break the degeneracy might be to hnd the Bayesian complexity of 
both models as in [16]. We expect such an analysis would be more punishing to the curvaton 
due to the additional parameters it requires. 

We hnd the required measurements and their associated precisions (t(/nl) to render 
the curvaton scenarios the same as the AGDM concordance model in the evidence, assuming 
that all other observational data remains the same. The conclusions from this analysis are 
twofold: hrstly, if /nl is as large as currently allowed at 2 — cr from Planck 2015, then Euclid 
should reach sufficient accuracy to favour the curvaton scenario and secondly, in the limit of 
the dominant curvaton scenario which predicts /nl = —5/4 an error bar signihcantly below 
1 will be required, which is forecast to be achievable in the longer term [21-24]. Unless /nl 
is signihcantly larger than unity, the error bar on /nl needs to be decreased by over an order 
of magnitude before non-Gaussianity has a signihcant effect on model selection. 

We have shown it is possible to hnd a simple multi-held model that is indistinguish¬ 
able from the single-held quadratic inhaton using the Bayesian evidence combined with the 
Planck 2015 data, despite having thee additional model parameters. The curvaton scenario 
is favoured when using the BIGEP/Keck/P/ancA: data. We suggest that this result may be 
viewed in two ways: hrstly, we may further strengthen the incentive to improve the current 
constraint on primordial non-Gaussianity as a tool to distinguish between inhationary poten¬ 
tials and secondly, one may take this result as motivation to explore the parameter space and 
Bayesian evidence of other multi-held models in future work. 

Acknowledgments 

RH was supported by a Royal Society Summer Research Studentship for the main duration 
of this work. GB is supported by a Royal Society University Research Fellowship. The 
authors wish to thank Tommaso Giannantonio, Antony Lewis, Andrew Liddle, Jerome Martin, 
Hiranya Peris, Ghristophe Ringeval and especially Vincent Vermin for helpful discussions. 

References 

[1] D. H. Lyth and D. Wands, Generating the curvature perturbation without an inflaton, 

Phys.Lett. B524 (2002) 5-14, [hep-ph/0110002]. 

[2] K. Enqvist and M. S. Sloth, Adiabatic CMB perturbations in pre - big bang string cosmology, 

Nucl.Phys. B626 (2002) 395-409, [hep-ph/0109214]. 


- 13 


[3] T. Moroi and T. Takahashi, Effects of cosmological moduli fields on cosmic microwave 
background, Phys.Lett. B522 (2001) 215-221, [hep-ph/0110096]. 

[4] S. Mollerach, Isocurvature Baryon Perturbations and Inflation, Phys.Rev. D42 (1990) 313-325. 

[5] A. D. Linde and V. F. Mukhanov, Nongaussian isocurvature perturbations from inflation, 
Phys.Rev. D56 (1997) 535-539, [astro-ph/9610219]. 

[6] K. Dimopoulos, D. Lyth, A. Notari, and A. Riotto, The Curvaton as a 
pseudoNambu-Goldstone boson, JHEP 0307 (2003) 053, [hep-ph/0304050]. 

[7] D. Langlois and F. Vernizzi, Mixed inflaton and curvaton perturbations, Phys.Rev. D70 (2004) 
063522, [astro-ph/0403258]. 

[8] T. Moroi, T. Takahashi, and Y. Toyoda, Relaxing constraints on inflation models with 
curvaton, Phys.Rev. D72 (2005) 023502, [hep-ph/0501007]. 

[9] T. Moroi and T. Takahashi, Implications of the curvaton on inflationary cosmology, Phys.Rev. 
D72 (2005) 023505, [astro-ph/0505339]. 

[10] K. Ichikawa, T. Suyama, T. Takahashi, and M. Yamaguchi, Non-Gaussianity, Spectral Index 
and Tensor Modes in Mixed Inflaton and Gurvaton Models, Phys.Rev. D78 (2008) 023513, 
[arXiv; 0802.4138]. 

[11] K. Dimopoulos, K. Kohri, D. H. Lyth, and T. Matsuda, The inflating curvaton, JGAP 1203 
(2012) 022, [arXiv: 1110.2951]. 

[12] N. Bartolo and A. R. Liddle, The Simplest curvaton model, Phys.Rev. D65 (2002) 121301, 
[astro-ph/0203076]. 

[13] C. T. Byrnes, M. Cortes, and A. R. Liddle, Gomprehensive analysis of the simplest curvaton 
model, Phys.Rev. D90 (2014) 023523, [arXiv; 1403.4591]. 

[14] Planck Collaboration Collaboration, P. Ade et ah, Planck 2015. XX. Gonstraints on 
inflation, arXiv; 1502.0211. 

[15] Planck Collaboration Collaboration, P. Ade et ah, Planck 2013 results. XXII. Gonstraints 
on inflation, arXiv; 1303.5082. 

[16] J. Martin, C. Ringeval, R. Trotta, and V. Vennin, The Best Inflationary Models After Planck, 
JGAP 1403 (2014) 039, [arXiv; 1312.3529]. 

[17] R. Easther and H. V. Peiris, Bayesian Analysis of Inflation II: Model Selection and Gonstraints 
on Reheating, Phys.Rev. D85 (2012) 103533, [corXiv; 1112.0326]. 

[18] T. Giannantonio and E. Komatsu, Bayesian evidence of nonstandard inflation: Isocurvature 
perturbations and running spectral index, Phys.Rev. D91 (2015), no. 2 023506, 

[arXiv; 1407.4291]. 

[19] G. Lazarides, R. R. de Austri, and R. Trotta, Gonstraints on a mixed inflaton and curvaton 
scenario for the generation of the curvature perturbation, Phys.Rev. D70 (2004) 123527, 
]hep-ph/0409335]. 

[20] J. Ellis, M. Fairbairn, and M. Sueiro, Rescuing Quadratic Inflation, JGAP 1402 (2014) 044, 
[arXiv; 1312.1353]. 

[21] T. Giannantonio, C. Porciani, J. Carron, A. Amara, and A. Pillepich, Gonstraining primordial 
non-Gaussianity with future galaxy surveys, Mon.Not.Roy.Astron.Soc. 422 (2012) 2854-2877, 

[arXiv; 1109.0958]. 

[22] S. Camera, M. G. Santos, and R. Maartens, Probing primordial non-Gaussianity with SKA 
galaxy redshift surveys: a fully relativistic analysis, cirXiv; 1409.8286. 

[23] O. Dore, J. Bock, P. Capak, R. de Putter, T. Eifler, et ah, SPHEREx: An All-Sky Spectral 
Survey, arXiv; 1412.4872. 


14 - 


[24] B. Leistedt, H. V. Peiris, and N. Roth, Constraints on Primordial Non-Gaussianity from 
800000 Photometric Quasars, Phys.Rev.Lett. 113 (2014), no. 22 221301, [arXiv; 1405.4315]. 

]25] D. Wands, N. Bartolo, S. Matarrese, and A. Riotto, An Observational test of two-field 
inflation, Phys.Rev. D66 (2002) 043520, [astro-ph/0205253]. 

]26] P. Creminelli and M. Zaldarriaga, Single field consistency relation for the 3-point function, 
JCAP 0410 (2004) 006, [astro-ph/0407059]. 

[27] J. Martin, H. Motohashi, and T. Suyama, Ultra Slow-Roll Inflation and the non-Gaussianity 
Consistency Relation, Phys.Rev. D87 (2013), no. 2 023514, [arXiv: 1211. 0083]. 

[28] X. Chen, H. Firouzjahi, M. H. Namjoo, and M. Sasaki, A Single Field Inflation Model with 
Large Local Non-Gaussianity, Europhys.Lett. 102 (2013) 59001, [arXiv: 1301. 5699]. 

[29] S. Mooij, G. A. Palma, and A. E. Romano, Consistently violating the non-Gaussian consistency 
relation, arXiv: 1502.0345. 

[30] A. D. Linde and V. Mukhanov, The curvaton web, JCAP 0604 (2006) 009, 

[astro -ph/051 1736] . 

[31] Planck Legacy Archive , base_ r_ WMAP, 2013. 

[32] R. Trotta, Bayes in the sky: Bayesian inference and model selection in cosmology, 
Contemp.Phys. 49 (2008) 71-104, [airXiv: 0803. 4089]. 

[33] S. Nesseris and J. Garcia-Bellido, Is the Jeffreys’ scale a reliable tool for Bayesian model 
comparison in cosmology?, JCAP 1308 (2013) 036, [arXiv: 1210.7652]. 

[34] H. Jeffreys, An Invariant Form for the Prior Probability in Estimation Problems, Proceedings 
of the Royal Society of London A: Mathematical, Physical and Engineering Sciences 186 
(1946), no. 1007 453-461. 

[35] K. Enqvist, R. N. Lerner, O. Taanila, and A. Tranberg, Spectator field dynamics in de Sitter 
and curvaton initial conditions, JCAP 1210 (2012) 052, [arXiv: 1205.5446]. 

[36] Planck Collaboration Collaboration, P. Ade et al., Planck 2015 results. XVII. Constraints 
on primordial non-Gaussianity, arXiv: 1502.0159. 

[37] D. Freedman and P. Diaconis, On the histogram as a density estimator:l 2 theory, Zeitschrift 
fAijr Wahrscheinlichkeitstheorie und Verwandte Gebiete 57 (1981), no. 4 453-476. 


- 15 


